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We investigate high-harmonic generation in closed sys- 
tems, using the two-level atom as a simplified model. By 
means of a windowed Fourier transform of the time-dependent 
dipole acceleration, we extract the main contributions to this 
process within a cycle of the driving field. We show that the 
patterns obtained can be understood by establishing a par- 
allel between the two-level atom and the three-step model. 
In both models, high-harmonic generation is a consequence 
of a three-step process, which involves either the continuum 
and the ground state, or the adiabatic states of the two-level 
Hamiltonian. The knowledge of this physical mechanism al- 
lows us to manipulate the adiabatic states, and consequently 
the harmonic spectra, by means of a bichromatic driving field. 
Furthermore, using scaling laws, we establish sharp criteria 
for the invariance of the physical quantities involved. Conse- 
quently, our results can be extended to a broader parameter 
range, as for instance those characteristic of solid-state sys- 
tems in strong fields. 



I. INTRODUCTION 

The generation of high-order harmonics of a strong 
laser field (/ ^ 10^"'W/cm^) in gaseous samples, where 
coherent light in the extreme ultraviolet regime is ob- 
tained from infrared input radiation, originated a break- 
through in nonlinear optics. In these systems, composed 
by atoms or small molecules, high-harmonic generation 
(HHG) is a well-understood issue |]. These highly non- 
linear spectra exhibit very particular features: a fre- 
quency region with harmonics of roughly the same in- 
tensities, the "plateau" , and a sharp decrease in the har- 
monic yield at the plateau's high-energy end, the "cut- 
off". Since the early nineties, not only these features 
have been investigated, but also the HHG time profile 
physical mechanisms and the propagation of 

the harmonic radiation in gaseous media . These stud- 
ies culminated with countless proposals of how to control 
high harmonics, as diverse as for instance polychromatic 
10-^ or static |lO| fields, ultrashort pulses or ad- 
ditional potentials Jl2| , many of them having even been 
realized experimentally |]l3| . 

One of the first models proposed to describe high- 
harmonic generation in atoms or diatomic molecules was 
a two-level atom Within this framework, a par- 

ticularly important paper is [p^ . Therein, it is shown 
that these harmonics are a consequence of the popula- 
tion transfer between the field-dependent states obtained 
from the diagonalization of the two-level Hamiltonian. 



This physical mechanism has not been investigated in de- 
tail, and there is a very simple reason for this apparent 
lack of interest: it turned out that an at first sight com- 
pletely different physical picture is far more successful in 
explaining high- harmonic generation for these systems. 
This picture, known as "the three-step model" , portraits 
high-harmonic generation as a process in which an elec- 
tron leaves an atom at an instant to (the first step) , prop- 
agates in the continuum being accelerated by the field 
(the second step), and recombines with the ground state 
of its parent ion |^ at a later time ti , emitting a high- 
harmonic photon (the third step). This model has shown 
that the interplay between a bound state and the contin- 
uum, which is not present in a two-level atom, is essential 
for a correct physical description of high-harmonic gener- 
ation. Thus, the three-step model has established itself as 
the paradigm for describing this phenomenon (see, e.g., 
for a comparison of both models) . 

Until very recently, only gaseous systems were believed 
to be possible high-harmonic sources, due to the high in- 
tensities involved. However, nowadays, this picture has 
changed. With the advent of short pulses, there are solid- 
state materials which can survive the necessary intensity 
regime, namely 10^^ — 10^'*W/cm^ This has led to 
theoretical studies on high-harmonic generation in ma- 
terials such as thin crystals |]l7| , or carbon nanotubes 
p8[ . Another example of a new and unexpected effect is 
for instance carrier-wave Rabi flopping, which has been 
recently measured experimentally p9| ] . 

Furthermore, apart from this entirely new parame- 
ter range, even for considerably lower driving- field in- 
tensities, as for instance / '-^ lO^W/cm^, one may in 
principle extend the frequency of far-infrared radiation 
{uj ~ IGHz) in up to two orders of magnitude by using 
adequate materials. For instance, for GaAs/A^Gai-ajAs 
wells intersubband transitions of loq ~ ITHz may serve 
this purpose [Q. Apart from these solid-state materials, 
HHG involving larger molecules is becoming a problem 
of interest |^,p2| . 

For these complex systems, it is not entirely clear 
whether bound-to-continuum transitions still yield the 
most adequate description of high-harmonic generation. 
In fact, recent studies have shown that, for aromatic 
molecules, transitions involving solely bound states are 
far more important for high-harmonic generation than 
the interplay between the ground state and the contin- 
uum p2[ . Thus, theoretical approaches in which the con- 
tinuum is not taken into account may be possibly used 
to describe this phenomenon in systems as, for instance. 
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quantum wells ]2C| , |23| -|26| . Furthermore, descriptions of 
nonlinear optical processes in solids are widely based on 
the Hartree-Fock semiconductor Bloch equations. Under 
special conditions, such as low doping density, equal effec- 
tive masses in both subbands involved, parallel subbands, 
and not too wide wells, these equations are formally iden- 
tical to those describing the evolution of a two-level atom. 
Otherwise, collective effects must be taken into account 
and this analogy is lost ||2^,|2|-||l. 

A common characteristic of all the above-stated sys- 
tems is their intrincated internal structure, with the pres- 
ence, as the external parameters are varied, of several 
level crossings. In particular concerning HHG, the peri- 
odic level crossings caused by the temporal dependence 
of the laser field are very important |14[| . Thus, in order 
to control the harmonic spectra also in this context, one 
needs to understand the interplay between the popula- 
tion transfer at these crossings and high-harmonic gen- 
eration. 

Even in the simplest case for which these level cross- 
ings occur, namely a two-level atom, it is only clear that 
most of the population transfer between the field-dressed 
states takes place at the level crossings. However, this 
does not necessarily mean that the population transfers, 
within a field cycle, which contribute to the generation 
of a particular group of harmonics occur at the level- 
crossing times. Unanswered questions in this framework 
concern not only these times, but also how they depend 
on the external-field parameters, such as its intensity and 
frequency, and how one can use this information to con- 
trol the emission spectra of a "closed" , non-ionizing sys- 
tem. Another interesting issue concerns the existence of a 
one-to-one correspondence between the three-step model 
and the two-level atom. This was proposed in due 
to the different time scales involved in the process, and 
in due to a formally identical expression describing 
population transfers in both models. In these references, 
however, there is no proof that this correspondence really 
holds. 

The answer to these questions is the main objective 
of this work. The paper is organized as follows: in Sec. 
H we briefly discuss the theoretical background for the 
studies performed in this paper. In the following sections 
we present our results. In Sec. [II, we concentrate on a 



detailed analysis of the population transfers and the time 
profile of harmonic generation for a monochromatic field. 
Subsequently (Sec. [V), we provide concrete examples of 
how an additional driving field may alter the periodic 
level crossings, and consequently the harmonic emission 
of a closed system. Furthermore, we address the scal- 
ing behavior of the physical quantities involved (Sec. ^), 
establishing sharp criteria for their invariance. Finally, 
in Sec. VI we close the paper with some concluding re- 
marks. 



II. BACKGROUND 



A. Two-level atom 



The simplest case for which level crossings occur, and 
a widely used approximation for describing physical sys- 
tems, is a two- level atom Within this picture, the 
time-dependent wave function is given by 



\m)^c^{t)\M + c,{t)\c^, 



(1) 



where C„(t) = ((/>„ \i'{t)) denotes the overlap of the total 
wave function with the n-th state of an arbitrary basis. 
The evolution of the system is described by the time- 
dependent Schrodinger equation, 



dt \ Ci(t) J - ^ [Ci{t) 



(2) 



where H is the Hamiltonian matrix, which, in our case, 
describes an atom in an external laser field. We use 
atomic units throughout. The basis states \cf>n) are cho- 
sen according to the problem at hand. We are par- 
ticularly interested in a basis which yields sharp, well- 
separated level crossings in the strong-field regime. 

A widely used basis are the field-free-states, also known 
as the "diabatic basis" . In this case, the Hamiltonian is 
given by 



-cjio/2 xiaE{t) 
xiQE{t) ujio/2 



(3) 



where ujiq is the transition frequency between the field- 
free bound states, E{t) = Eof{t) is the external field 
and xio the dipole matrix element ^(/)^| x ), where 
\(f)n) denotes the field- free, "diabatic" basis states. This 
basis is very convenient for studying level crossings in the 
low-intensity laser field regime. For strong laser fields, 
however, the field-free states are too strongly mixed, such 
that a more appropriate basis is needed. Such a basis, 
which will be called by us "exchanged basis" , is obtained 
applying the unitary transformation 



Ud- 
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1 1 

-1 1 



(4) 



onto the diabatic basis. The transformation (Q) was 
used in to interchange the diagonal and the non- 
diagonal terms of the Hamiltonian (0). In this case, the 
exchanged-basis energies = ±xioE{t) cross, and the 
coupling which causes the crossing is effectively given 
by ljiq/2. The crossings occur within a time interval 
to — tc < t < to -\- tc-, where tc is the time for which 
the off-diagonal and diagonal terms of the Hamiltonian 
become equal and is the time for which the field van- 
ishes. For strong enough fields, the times over which the 
crossings take place are much smaller than the period of 
the driving field. Thus, to first approximation, one may 
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assume that the crossings take place instantaneously at 
to. In the following we call to "crossing times". 

Another important set of basis states are these which 
diagonalize H. This basis is the so-called "adiabatic ba- 
sis", in the sense that the states "follow" the field, and 
is obtained by means of the unitary transformation 

Un^A=(''''^ "'"^V (5) 
\ — sm X cos X J 

with X = — l/2arctan(22;io£'(i)/ijJio)- This gives 

if^ = Ud^aHUE^a = ( 0^ °A ) , (6) 
where the field-dressed energies are given by 

si=±yu;l, + i2x^oEit))^. (7) 

Applying U o-tA to the diabatic basis states, one obtains 
the field-dressed, "adiabatic" states 

=cosx|<^?)+sinx|0f) (8) 

and 

|</)f(0) = -sinx|0?)+cosx|<^f), (9) 

whose energies are, respectively, and . In order 

to compute the harmonic spectra, one needs the Fourier 
transform of the time-dependent dipole. This quantity is 
given, in its length and acceleration form, by 

x = xio[g{t)cos2x + h{t)sm2x], (10) 

and 

X — —ij^l^x + 2LUiQxlQE{t) [h{t) cos 2x — g{t) sin 2x\ , 

(11) 

respectively, with g{t) = C^'^{t)Cf{t)+Cl'^{t)C^{t) and 
h{t) = \C^{t)\' - \C^{t)\', where C^{t) = {c^^{t) \m) 
denotes the projection of the wave function |V'(t)) onto an 
adiabatic state. The equations above are the superposi- 
tion of two distinct terms, namely the crossed terms and 
the population difference between the adiabatic states. 
Since the population difference h{t) roughly "follows" the 
field, it contributes mainly to the generation of low har- 
monics, whereas g(t) is expected to be responsible for the 
high harmonics. This has been confirmed by numerical 
studies (not shown). 

An interesting feature is that, in the extreme limit 
Eq —> oo, the transformation (H) formally corresponds 
to (||) and the dipole length ( p^ becomes proportional 
to the population difference between the adiabatic states. 
However, one should keep in mind that, only in this limit, 
the states obtained using (||) on the field-free states and 
the adiabatic states are formally equivalent. In general, 
this is not the case. 



In the subsequent sections, we work mainly in the adia- 
batic basis, and refer to crossings of the exchanged-basis 
energies. For the adiabatic energies, there are avoided 
crossings. The results discussed in this paper have been 
obtained from the numerical solution of Eq. (^) in the 
adiabatic basis, by means of a fourth-order Runge-Kutta 
method. Unless stated otherwise, the driving field is 
turned on instantaneously. 

B. Windowed Fourier transform 

For both open and closed systems, high-harmonic gen- 
eration is always related to abrupt population transfers. 
Depending on the group of harmonics, they occur at 
particular times, which give the main contributions to 
high-harmonic generation within a field cycle. For an 
atom in a strong laser field, for instance, these times are 
well-known and correspond to the return times ti of an 
electron which left an atom at a previous time Iq. For 
a closed system, the times to correspond to the level- 
crossing times and the times ti are still an open question 
to some extent. A very useful method to extract these 
latter times from the time-dependent dipole is performing 
a Fourier transform with a temporally restricted window 
function. For an arbitrary function f(t'), this transform 
is 

-t-oo 

^(t,0,(T)= J dt'J{t')W{t,t\n,a) , (12) 

— oo 

where i, U, and a denote the time and harmonic frequency 
at which the window function is centered, and its tempo- 
ral width, respectively. We consider a Gabor transform, 
for which the window function is given by 

W{t, t', n, cr) cxp[-(t - t'f/(7^] exp[im'] . (13) 

The usual Fourier transform which yields no tem- 

poral information, is recovered for cr ^ oo. The temporal 
width a corresponds to a frequency bandwidth ctq = 2/ a. 
For a temporal width smaller than the period T — 211/0; 
of the driving field, the peaks in the time-resolved spec- 
tra |jF(i, rj, cr)p yield the recombination times ti. This 
method has been extensively used in the literature, in 
the three-step model framework 

III. GENERAL PICTURE 

We shall now investigate the connection between HHG 
and the periodic level crossings in detail and draw a gen- 
eral physical picture of the mechanisms involved. The 
simplest physical situation for which one can do this is a 
monochromatic field 

E{t) = Eosin{ujt), (14) 
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where Eq and ui denote the field strength and frequency, 
respectively. In this case, the time tc is given by the 
condition 



2x10^^0 ' 



(15) 



If the field amplitude Eq is large enough, then ojtc <C 1, 
and the avoided crossings of the adiabatic states are well- 
separated. Thus, the crossing times to are well-defined 
and there is efhcient population transfers at to- Hence, 
one expects the corresponding spectrum to exhibit a wide 
plateau and a sharp cutoff. 

The avoided crossings occur at the times tg = m:/uj 
for which the field is vanishing. Thus, one expects the 
population transfers between the states \4'n{t)) to occur 
at these times. This is partially confirmed by Fig. 1, 
where the populations of the adiabatic states are plotted 
as functions of time. In fact, the pronounced peaks at 
the times Iq clearly show that most population transfer 
takes place at these times. There are however several 
smaller peaks, which are symmetric with respect to the 
times tiM = i'2n + l)7r/2w for which the field is maximal. 
These peaks show that population transfer also occurs at 
other times, and can be seen in detail in Fig. 1(b). 
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FIG. 1. Populations |C^(t)p of the adiabatic states as 
functions of time, for transition frequency toio ~ 0.409 a.u., 
external field parameters uj = 0.05a.u., Eq = 0.6 a.u., and 
dipole-matrix element a;io = 1.066 a.u.. Part (a) shows this 
feature for several cycles of the driving field, whereas part 
(b) depicts the population of the excited adiabatic state only 
within half a cycle. The times are given in units of the field 
cycle T = 2n/Lj. The driving field is turned on linearly within 
two periods. 



The role of these population transfers in HHG can be 
understood using the Gabor transform of the dipole ac- 
celeration. The peaks in the Gabor spectra give the main 
contributions for high-harmonic generation within a field 
cycle. For the cutoff harmonic, there is a single peak 
at tiM which splits into two, for the plateau harmonics. 
This peak gets further apart as the harmonic frequency 
decreases, varying from tiM to the times at the imme- 
diate vicinity of the avoided crossings. These results are 
displayed in Fig. 2. 




6.00 



FIG. 2. Gabor spectra of the dipole acceleration (Eq. (11)) 
as functions of time, for field strength Eq — 1 a.u., field fre- 
quency Lu = 0.05 a.u., transition frequency ti;io — 0.409 a.u. 
and dipole matrix element xio = 1.066 a.u.. The cutoff har- 
monic lies at Qm ~ 43cj. The time width of the window func- 
tion was chosen a — O.IT. Its center was chosen at the cut- 
off harmonics, as well as at harmonic energies which roughly 
correspond to = O.SJIm, = O.GJIm, and = OAi^M- AH 
time-resolved spectra have been normalized. The times are 
given in units of the field cycle T — 2tt/uj. The driving field 
is turned on linearly within two periods. 

The physical interpretation of these features is rather 
simple. At the times the level crossings occur, i.e., at to = 
nT/2, there is a population transfer from the adiabatic 
state |</'o (t)) to l^i^(t)) . The system remains in \(f>i{t)) 
until a further time ti, decaying back to |0^(i)) and 
emitting a harmonic of frequency = Nuj — e^— e^. The 
explicit expression relating the time ti to the harmonic 
frequency would then be 



Luti = arcsin ± V(A^7i)^ - (72)^ 



(16) 



with 71 = ti)/(2a;io£'o) and 72 = a;io/(2xio£'o). The 
physical significance of_7i and 72 will be discussed later 
in this paper (Sec. In order to obtain a har- 

monic at the maximum possible frequency Qm (i-e., 
the cutoff harmonic), the population transfer between 
the time-dependent states must occur at the times for 
which the energy difference — is maximal, i.e., at 
tiM = (2n -I- l)7r/2u;. As the harmonic energy decreases, 
there are two possible times for this population transfer 
to occur, a shorter and a longer one. The interference 
between these two possible quantum paths originates the 
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two-level atom plateau. This process repeats itself every 
half cycle of the driving field. This picture is supported 
by the fact that all peaks in the time-resolved spectra 
satisfy Equation (jl^) and thus can be traced back to 
population transfers between the adiabatic states. The 
times given by for the parameters of Fig. 2, together 
with the corresponding harmonic energies, are written in 
Table 1. 

An analogous picture is observed within the three-step 
model framework. The cutoff harmonic can only be gen- 
erated by an electron which returns to its parent ion 
with maximal kinetic energy. This maximal energy cor- 
responds to a particular return time, which appears as 
a single peak in the Gabor yield. Within the plateau, 
there are two possible sets of electron trajectories cor- 
responding to the same harmonic energy, such that this 
single peak splits into two In our case, the "first step" 
would be the population transfer from \(j)o{t)) to \4>i'{t)') 
at ig, the "second step" would be the system following 
adiabatically in a time interval t ~ ti ~ to and 
the "third step" the population transfer from [(/'^(i)) to 
\<j>o{t)') at ti, with subsequent harmonic generation. The 
corresponding physical picture is illustrated in Fig. 3. 
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FIG. 3. Schematic representation of high-harmonic genera- 
tion in a two-level atom. The population transfers at the level 
crossings occur at the times to and the main contributions to 
HHG occur at the times ti. The times tiM, tn and ti2 corre- 
spond to the generation of the cutoff and plateau harmonics, 
respectively. The main physical processes are indicated by 
arrows in the figure, and the corresponding energies can be 
read in the vertical axis. The adiabatic energies are given in 
units of the maximal energy e^j and the time in units of the 
field cycle. The field parameters are chosen in such a way 
that the ratio between the cutoff energy Qm — 2£^j and the 
transition frequency is 57m/wio = 10. 

Another interesting feature is that the population 
transfers between the adiabatic states are not strictly pe- 
riodic within t:/lu. Indeed, superposed to them, there are 
oscillations which occur within much larger time scales, 
their periods comprising several cycles of the driving field 
pof . These oscillations are also present in the dipole 
length and acceleration as a gobal enveloping function, 
whose amplitude, form and periodicity depend on the 



field strength Eq, the field frequency lu and on the dipole 
matrix element Xiq in a non-trivial way. These structures 
seem not to influence the harmonics gobally, but mainly 
the substructure of the spectra and the hyper-Raman 
lines (2|. 

In Fig. 4, we show these enveloping functions for the 
populations of the adiabatic states (Fig. 4(a)), the dipole 
acceleration (Fig. 4(b)), and the Gabor spectra of the 
plateau and cutoff harmonics (Fig. 4(c)). One should 
note that this enveloping function is the same for the 
Gabor transforms of all groups of harmonics displayed. 
Furthermore, it does not affect the splitting of the peaks, 
such that the population transfer times are always given 
by Eq. (©. 
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FIG. 4. Global structures as functions of time, for: (a) the 
populations |C^(t)p of the adiabatic states; (b) the dipole 
acceleration x{t); (c) the Gabor spectra of the cutofi^ and 
plateau harmonics. The time width of the window func- 
tion is a — O.IT. The field strength, the field frequency, 
the transition frequency and the dipole matrix element were 
chosen as Eq = 0.6 a.u., lj = 0.05 a.u., wio = 0.409 
a.u. and xio ~ 1.066 a.u., respectively. These parameters 
give 7i — 0.0391, 72 = 0.3197 and a cutoff frequency at 
f^max ~ 27uj. All Gabor spectra have been normalized to 
the maximum value obtained with the window function cen- 
tered at the cutoff. The field is turned on linearly within two 
periods. The time is given in units of the field cycle. 
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IV. CONTROL 

In this section, we consider a bichromatic driving field 

E(t) = Eoi sin(wt) + E02 sm{nujt + 6), (17) 

with two main purposes. First, we wish to confirm the 
physical picture in which the main contributions to a 
particular set of harmonics, within a field cycle, occur at 
the times ti such that the corresponding harmonic fre- 
quency is the difference — between the energies of 
the adiabatic states. Second, we are interested in under- 
standing how an additional field can be used to distort 
the avoided crossings between the adiabatic states in such 
a way that the harmonic emission can be controlled. In 
the bichromatic case, depending on the field parameters, 
the spectra may have several cutoffs, which are given by 
the maxima of — e;^. Consequently, the main contribu- 
tions to the generation of the cutoff harmonics take place 
at the times tiM for which these maxima occur. 

In order to obtain the level-crossing times to, as well as 
the times tiM, one needs the extrema ±sfj of the field- 
dressed energies e^. For the bichromatic field (|l^) they 
are given by 

cos(wi) -I- nC cos(nwt + 9) = (18) 

and 

sin{ujt) + C sm{nujt -f 6*) = 0, (19) 

where ( — E02/EQ1 denotes the field-strength ratio. 
Equation ( [l8| ) gives the extrema which coincide with 
those of the field, and therefore iiM, whereas Equation 
( |l9| ) gives those which correspond to the avoided cross- 
ings, and therefore Iq. Depending on the frequency ratio 
n, the field-strength ratio ( and the relative phase 6, these 
times, as well as the corresponding extrema, can be very 
different. In this paper, we will provide concrete exam- 
ples for a w — 2a; field, i.e., with n = 2, relative phases 
di = and 02 — JjJ'^j and arbitrary (. For these spe- 
cific parameters, (|l8|) and (|l^) have a simple form, with 
analytical solutions. 

A. 9^0 

In this case, Eq. ( |l8|) reduces to 

cos'^ {cut) + ^ cos{ujt) - ^ = 0, (20) 

which yields two sets of times, namely 

trM = ;^arccos(-l±iy^I^). (21) 

The solutions corresponding to the positive root exist for 
all field-strength ratios, whereas the remaining solutions 



are only present for C > 0.5. Further in this section, it 
will be shown that the first set gives the absolute maxima 
of e^, which correspond to the cutoff in the harmonic 
spectra, whereas the second set yields local maxima at 
much lower energies. 

The expression giving the avoided crossings, on its 
turn, can be written as 

sin(wi) [1 + 2C cos(wt)] = 0. (22) 

This equation yields the crossing times to = nTr/co, and 
= l/ui arccos[— 1/(2^)]. The crossing times to do not 
depend on the field-strength ratio and are the same as 
in the monochromatic case, whereas the crossing times 
t'o clearly do. Furthermore, these latter times are only 
present for ( > 0.5. 
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FIG. 5. Energies of the adiabatic states for a bichromatic 
field E{t) = £01 sm{ujt) + E02 sm{2ujt + 9), for 6 ^0 and sev- 
eral field-strength ratios ( — E02 / Eoi . The time t is given in 
units of the field cycle T = 2tt/u) and the field-dressed ener- 
gies in units of the maximal energy Em. The field parameters 
were chosen such that ^Im /uio = 8. 

Fig. H gives concrete examples of how the adiabatic 
energies depend on time, for different field-strength 
ratios. In contrast to the monochromatic case, is not 
periodic within half a cycle of the driving field. This is 
not surprising, since the periodicity of the field-dressed 
energies is effectively determined by E'^{t) (c.f. Eq. (^). 
For a monochromatic field, E'^{t) = E'^{t + tt/uj) always 
holds, whereas in the bichromatic case this is only true 
for odd frequency ratios n. This is clearly not the case 
addressed in this paper. For the phase = 0, one ob- 
serves that £±(t) = e^{2-K I uj ~ t) , if both times are taken 
symmetrically with respect to tg = mi/uj. This property 
already reflects itself in the expressions for toj tiM and 

derived in this section. 

Furthermore, one clearly sees that, as predicted in Eq. 
(p^), for C, < 0.5, the second driving wave only distorts 
the avoided crossings, making them broader at — (2n-|- 
l)7r/a; and sharper at to = 2m: /lo. For C, = 0.5, the broad 
crossing starts to split, originating the crossings given at 
the times t'^. This splitting also leads to the second set 
of maxima predicted by Eq. (pl|), which corresponds to 
a set of harmonics of relatively low frequencies. 
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One must now understand which consequences this ef- 
fect has on the physical quantities involved. With that 
purpose, we choose the strengths of both driving waves 
such that efj, and therefore the cutoff energy, remains 
unchanged and is equal to the monochromatic cutoff en- 
ergy, for variable field-strength ratio This gives 



Eoi — 



with /3 — cos(tiM)- 



En 



71^(1 + 2/30^ 



(23) 




FIG. 6. Populations |C^(t)|^ of the adiabatic states (Part 
(a)) and dipole acceleration (Part (b)) as functions of time, 
for a bichromatic field E{t) — Eoi sm{u!t) + E02 sin(2ajt -I- 9), 
with e = 0, LJ = 0.05 a.u., ujw = 0.409 a.u., xio = 1.066 
a.u. and field-strength ratio — Eq^/Eqi = 0.5. The field 



amplitudes were chosen according to Eq. (|23[), with Eo — 1 
a.u.. The time t is given in units of the field cycle. 

The population transfers between the adiabatic states, 
as functions of time, also exhibit very similar asymme- 
tries to the ones observed in the field-dressed energies. 
The population transfers at the broad crossings, for in- 
stance, take place at longer time intervals than those at 
the sharp crossings, making the oscillations in |C^(t)p 
asymmetric with respect to the times tiM- This asymme- 
try increases with increasing C. An example is provided 
in Fig. ^(a) . A similar feature occurs for the dipole accel- 
eration. This highly oscillating function exhibits nodes 
at the level-crossing times. In the monochromatic case, 
these nodes extend over identical temporal regions ev- 
ery half-cycle of the driving field. For bichromatic fields, 
however, with the distortion of the crossings by the sec- 
ond driving wave, this picture changes. There exist nar- 
rower and broader nodal regions, corresponding to the 



narrower and broader crossings, respectively. Thus, the 
oscillations of the dipole acceleration get "squeezed" be- 
tween the broader nodes. This feature can be seen in Fig. 
1(b). 
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FIG. 7. Gabor spectra of the dipole acceleration as func- 
tions of time, for a bichromatic field E{t) — Eqi sm{ujt) 
-hSo2 sin(2tjt + 9), with 9 ^ Q, uj ^ 0.05 a.u., low = 0.409 
a.u., 2:10 = 1.066 a.u. and several field-strength ratios 
C = E02/E01. The maximal field strength is kept fixed ac- 
cording to Eq. (jisj), with Eq = 1 a.u. The cutofi' energy lies 
at 0,M = 2e^j = 43lj. The temporal width of the window 
function is cr = O.IT. In Part (a), the window function is 
centered at the cutoff harmonics, and the field-strength ratio 
is < < 0.8. In Part (b), the center of the window function 
is taken for different frequencies, and ( — 0.8. All curves in 
the figure have been normalized to their maximum values. 

The Gabor transform of the dipole acceleration, taken 
at the cutoff and in the plateau, confirms this picture. In 
Fig. 0(a), there is a clear displacement of the peaks in 
the time-resolved spectra for the cutoff harmonics, with 
respect to the monochromatic case, and these peaks oc- 
cur at the times predicted by Eq. (^). Similarly to the 
monochromatic case, these peaks split into two in the 
plateau region, being, however, slightly asymmetric (Fig. 
0(b)). This asymmetry is related to the above-mentioned 
difference in the shapes of the crossings. Furthermore, for 
larger field-strength ratio, the additional times can also 
be seen, for a group of harmonics at the low-energy end 
of the plateau. The times to and tiA/, together with the 
respective cutoff energies, are given in Table II for the 
specific parameters considered in this figure. 
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B. 9 = tt/2 

For this relative phase, Equation (|l^) has the form 

cos{Lut) [1 - 2C sin(wt)] = 0. (24) 

This equation has two types of solutions: tiM = {n + 
1/2)tt/uj, which do not depend on the field-strength ra- 
tio and yield the same maxima as in the monochromatic 
case, and t[j^,j — l/tjarcsin[l/(4C)], which clearly depend 
on ( and exist only for ( > 0.25. This already hints at a 
completely different situation as in the previous section, 
which will now be discussed in detail. This also holds 
for the times at which the avoided crossings occur. They 
must now satisfy 



1 



1 



sin (ujt) - — sin(wt) - - = 



such that 



1 



tf) = — arcsin — ± — 



OJ 



1 



1 



4C 



1 

4(2 



(25) 



(26) 



all of them depending on (. This means that, in contrast 
to the case 9 — 0, one may shift all level-crossing times 
by changing the relative intensities of the driving waves. 
The set of crossings given by the positive root in ( |2^ ) 
exists only for C ^ li whereas the remaining crossings 
occur for all C- 

In Fig. ^ we depict the adiabatic states as functions 
of time, for several values of C, similarly to what was 
done for 6 = 0. This figure illustrates how the relative 
phase can radically alter the whole physical picture. For 
9 = 7r/2, already a relatively weak high-frequency wave 
considerably distorts the avoided level crossings, as well 
as the maxima of the field-dressed energies. An interest- 
ing feature is that the avoided crossings now move with 
the field-strength ratio. Furthermore, the maximal ener- 
gies are no longer equal, but, within a field cycle, there 
are two comparable and different cutoff energies. This 
can be directly seen by computing the extrema of the 
energies e^, which occur for tiM- 

For field-strength ratio ( < 0.25, they give the energies 



and 



'Ah 



'10 



cjtn + 4a;2p(£'oi 



(27) 



(28) 



which correspond to the times ^imi = 0.25T mod T, and 
to tiM2 = 0.75r mod T, respectively. These times define 
symmetry axes for the time-dependence of the adiabatic 
energies. 

For C > 0.25, a further splitting of the set of maxima 
at tiMi occurs, as predicted in Eq. (p4|). There exist now 
two sets of maxima, at the times t[]^.j, whose energies are 
equal and given by 




4x2 (1 + 8C2)2 



(29) 



These maxima are symmetric with respect to ^imi • For 
these times, the adiabatic energies now exhibit a mini- 
mum. This causes, for large additional avoided cross- 
ings (c.f. Fig. 8 for C = 0.8). The population transfers 
at these times are however small, and play only a sec- 
ondary role in the problem addressed in this paper. For 
the sake of simplicity, even after the second splitting, we 
shall refer to the lower-energy set of maxima as ejlj_^ . The 
other set of maxima does not split, and the correspond- 
ing times tiM2 remain constant for all C,. One should note 
that the adiabatic energies, in the 9 = 7r/2 case, satisfy 
e^(i) = £^(T/2 ~ t), if both times are chosen symmet- 
rically with respect to tiM-i or tiM2 ■ This also holds for 
the population-transfer times derived in this section. 
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FIG. 8. Energies of the adiabatic states for a bichromatic 
field E{t) = Eoi sm{ujt) + E02 am{2ujt -I- £»), for 6* = 7r/2 and 
several field-strength ratios — Eo2/Eoi. The time t is given 
in units of the field cycle and the field-dressed energies in 
units of the maximal energy e^j^ . The field parameters were 
chosen such that 57^/2 /'^lo = 8. The times tiMi are indicated 
in the figure by the dotted and solid grid lines, respectively. 

In order to investigate how the distortions in the 
adiabatic-state energies influence the physical quantities 
of interest, we shall keep the cutoff energy Qm^ = "^^th 
fixed, and equal to the cutoff energy of the monochro- 
matic case. Thus, the field strengths Eqi and Eq are 
related by 



£^01 — 



En 



1 + C 



(30) 



As in the previous section, we can trace all distortions 
observed in these physical quantities back to those ob- 
served in time dependence of e^. For instance, the shifts 
in the level-crossing times to predicted by Eq. ( p6| ) are 
also present in the main population-transfer times for 
the adiabatic states (Fig. J9(a)) and in the nodes of the 
dipole acceleration (Fig. H(b)). Another effect which is 
clearly seen in both quantities is the splitting of the max- 
ima near tiMi = 0.25T mod T. Indeed, there exist now 
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two sets of maxima which are symmetric with respect to 
these times, for ( > 0.25. 




FIG. 9. Populations |C^(t)p of the adiabatic states (Part 
(a)) and dipole acceleration (Part (b)) as functions of time, 
for a bichromatic field E{t) — Eoi sm{uit) + E02 sin{2Ljt + 6), 
with e = 7r/2, uj = 0.05 a.u., wio = 0.409 a.u., xio = 1-066 
a.u., and field-strength ratio ( = E02/E01 = 0.8. The time t 
is given in units of the field cycle. 

We now investigate the Gabor transform of the cutoff 
and plateau harmonics. In Fig. 0(a), we display the 
time-resolved spectra, centered at the harmonic frequen- 
cies 17 Afa = 2e;^^^, for different field-strength ratios C- The 
monochromatic case is also displayed for comparison. As 
a general feature, for C 7^ 0, the peaks of the Gabor spec- 
tra at iiMi = 0.25T mod T vanish. This is a direct con- 
sequence of the splitting of the extrema of the adiabatic 
energies caused by the high-frequency wave. Due to this 
splitting, the energy maxima near tiMi lie outside the 
range of the window function and do not contribute to 
the time-resolved spectra. Furthermore, as predicted in 
Eq. (25), the peaks at the maxima iij\/2 = 0.75T mod T 
do not move in time as ^ is varied. 

Taking now the window function ( |l3|) centered at 
fi^/j — 2ef,j^ (Fig. ^(b)), one observes, as expected, a 
completely different behavior for the peaks near tiMi — 
0.25T mod T. For ( < 0.25, these peaks are exactly at 
these times. For ( > 0.25, as expected, they now occur 
at t[j^.j — arcsin[l/(4C)], which vary with the field- 
strength ratio (. Furthermore, this second set of peaks 
splits for these larger field-strength ratios, such that two 
sets of peaks which are symmetric with respect to tiMi 



are now present. Other sets of peaks which can be seen 
in the picture correspond to the upper-plateau return 
times, which occur for f2 < fl]\i2 and are symmetric with 
respect to tiM2 — 0.75T mod T. These peaks come from 
the splitting of tiM2 j which occurs in this energy range 
(c.f. Fig. H). The population-transfer times for the spe- 
cific parameters of this figure, together with the corre- 
sponding harmonic frequencies, are given in Table III. 
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FIG. 10. Gabor spectra of the dipole acceleration as a 
function of time, for a bichromatic field E{t) — Eoi sm{ujt) 
+E02 sm{2Lut + 6), with 9 ^ -it /2, cj = 0.05 a.u., ujio = 0.409 
a.u., xio — 1.066 a.u. and several field-strength ratios 
C = £■02 /Eoi . The maximal field strength is kept fixed accord- 
ing to Eq. ( po[ ) and equal to _Bo = 1 a.u. The upper-cutoff 
energy lies at ^1/2 = iSuj. The lower-cutoff energy varies 
with (. All cutoff energies are given in Table III, together 
with the population transfer times to and tiM^. In Part (a), 
the window function is centered at the upper-cutoff harmon- 
ics {Qm2 ~ 2eff2), and the field-strength ratio is < < 0.8. 
In Part (b), the center of the window function is taken at 
flM, = 2eij^. All curves have been normalized to their max- 
imum values. In Part (a), the monochromatic case is also 
displayed for comparison. 



C. Fourier spectra for the two phases 

In the investigations performed so far, our main objec- 
tive was to understand how an additional driving wave 
may distort the time dependence of the adiabatic en- 
ergies and the time profile of harmonic generation. In 
this section, we address the question of how these distor- 
tions influence the harmonic spectra. Furthermore, we 
are interested in extending the cutoff, and, by doing so, 
guaranteeing that the harmonics in this energy region 
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are strong enough for applicational purposes. Clearly, 
the ideal scenario is to extend the cutoff energy without 
any intensity loss in the corresponding harmonic range. 

With that purpose, we keep Eqi and Eq2 fixed and 
compare spectra obtained for 9i = and 62 — tt/2. These 
results are displayed in Fig. ^ As a global feature, one 
observes that, for 6* = 0, all harmonics behave in a very 
similar way, with no distinct regions, as for instance a 
double plateau, in the spectra. This is related to the 
fact that no splitting of the cutoff energy occurs in this 
case. The two maxima in have the same energy, even 
though the level-crossing pattern is no longer periodic 
in T/2. On the other hand, for 9 — tt/2, there is a 
clear double-plateau structure. In fact, one can identify a 
completely different physical behavior for the harmonics 
in the frequency regions fl < Qmi and ^Imi < ^ < ■ 
The double-plateau structure is due to the different cutoff 
energies which exist in the 6 = tt/2 case. 

Another generic feature is that the cutoff energy is ex- 
tended for = 7r/2. This is expected, since this quantity 
is given by the maximum energy difference between the 
adiabatic states. For a field given by Eq. (p7|), the maxi- 
mal possible energy is obtained for E{tiM2j = Eqi +^'02- 
This yields the harmonic frequency M2 -i discussed in the 
previous subsection. 

There exist however non-generic features, which de- 
pend on the absolute field parameters, as for instance its 
strength. Examples of such features are the intensity ra- 
tio between the upper and lower parts of the plateau for 
= tt/2, and the intensities of the harmonics obtained 
for 6 — tt/2, compared to those obtained for = 0. Thus, 
depending on the absolute parameters used, it is not al- 
ways possible to extend the cutoff energy without loss of 
intensity. In order to control HHG in a two-level atom, 
a more detailed study of these features for the particular 
system in question is necessary. 
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FIG. 11. Spectra computed from the dipole acceleration, 
for the bichromatic field E{t) — Eoi sin(tjt) +£02 sm{2ujt + 6), 
for 6 — 0, 6 = tt/2, and field strenghts -Eoi ~ 1-0 a.u. and 
£-02 ~ 0.2 a.u. The field is switched on linearly within two cy- 
cles. The remaining parameters are lu = 0.05 a.u., cjio = 0.409 
a.u., xio = 1.066. The cutoff frequency for 6 = is roughly 
at Qm — 46(jj, whereas for 6 = tt/2 the cutoff frequencies are 
approximately at Omi ~ 35aj and ^Im2 ~ 52uj. All cutoff 
energies are indicated by arrows in the figure. 

These non-generic features are mainly due to the fact 
that the population transfer at the level crossings is, in 
general, given by more complicated expressions than in 
the monochromatic case. Indeed, these expressions de- 
pend on the shape and width of the crossing, and on 
the duration of the interaction. These shapes have been 
studied in 32|. Furthermore, the global structures of 
the adiabatic-state populations |C^(i)p have a stronger 
influence on the spectra in the bichromatic case than for 
monochromatic driving fields. 



V. SCALING BEHAVIOR 

In the results discussed in the previous sections, we 
have used rather unrealistic frequencies and intensities 
for the driving fields, for which most physical systems 
would ionize immediately. This choice of parameters al- 
lows us to obtain results with very little numerical effort. 
In order to extend our computations to more realistic 
cases, as for instance solids, there are two possibilities. 
Either one slightly increases the effort to obtain the nec- 
essary precision, or one must find specific combinations 
of parameters for which the physical quantities involved 
remain invariant. This second approach has the advan- 
tage to provide additional insight into the physics of the 
problem. 

With that purpose, we analyze the scaling behavior 
of these quantities. We use scaling laws which have 
been derived elsewhere [Q, in the context of stabiliza- 
tion of atoms in strong laser fields. We concentrate on 
the question of whether driving fields of much lower fre- 
quencies and intensities could originate similar spectra, 
with, for instance, the same number of harmonics, or the 
same population-transfer times, in units of the field cycle. 
Therefore, our starting point will be the expression 



sin(a;ti) -I- (sm{nLjti + 9) = ± 



(72)' (31) 



which relates the harmonic energy to the energy differ- 
ence of the adiabatic states. This equation gives the 
population-transfer times. For C = 0, one has the 
monochromatic- field case (Eq. (f6h), and, for ( ^ 
and n = 2, the bichromatic situation discussed in the 
previous section. Note that the parameters Eq, lu, luiq 
and xiQ appear combined, as 71 = uj/(2xioEo), or 72 = 
luio/(2xiqEo). The denominators of these expressions 
give the Rabi frequencies — 2xioEq, which scale like 
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the energies (c.f. Eqs. (||) and (g) for the two-level Hamil- 
tonian). This keeps the Schrodinger Equation invariant 
under scale transformations. 

We now consider the scale transformation 



U! 



Lu' = Xuj; uJiQ 



'10 



n'„ = xn 



R, 



(32) 



where A denotes the dilatation factor. The invariance 
of the Schrodinger equation also requires that the time 
scales as t t' — X^^t, such that Eq. (|3|) will remain 
invariant. 

This apparently trivial result has far-reaching conse- 
quences. In fact, it shows that, for any set Eq, lu, wiq and 
xio, the number of harmonics N in the spectra and the 
corresponding population-transfer times ti = ujti/{27r), 
given in terms of field cycles, remain invariant, as long 
as 7i and 72 are kept constant. 

Since the unitary transformation which gives the 
adiabatic states also depends on Eq, uj, luiq and xio 
through 71 and 72, it also remains invariant in this case. 
Thus, this invariance must also hold for the populations 
of these states, i.e., \C^{t)\^ = \C^{t')\^. 

Another quantity of interest is the dipole acceleration. 
A quick inspection of Eq. (|ri ) shows that this quantity 
does not remain invariant under the above-stated trans- 
formations. In fact, it scales as xio multiplied by the 
square of the energy. The dipole matrix element scales 



as xio 



= Thus, xit) = X^/^x{t'). 




FIG. 12. Global structures as functions of time, for: (a) the 
populations |C^(t)p of the adiabatic states; (b) the dipole 
acceleration x{t). The field strength, the field frequency, 
the transition frequency and the dipole matrix element were 
chosen as Eq = 6.71 x 10~® a.u., tj = 2.5 x 10"^ a.u., 
cjio — 2.045 X 10~^ a.u., and a;io = 47.673 a.u., respectively. 
These parameters are typical for solid-state systems and give 
71 — 0.0391, 72 = 0.3197, which are the same as in Fig. 4. 
They are obtained from those in Fig. 4 using a scaling trans- 
formation with X — 1/2000. For this set of parameters, we 
have used a five times smaller timestep than in the previous 
figures and double precision. The dipole acceleration is given 
in atomic units and the time is given in units of the field cycle. 
The field is switched on linearly within two cycles. 

The above-stated conclusions are confirmed by Fig. |l^. 
In this figure, we display the same physical quantities as 
in Fig. 4 for a completely different set of parameters 
which, however, yield the same 71 and 79. The popu- 
lations \C^{t)\'^, in this case (c.f. Fig. |l^(a)) are, as 
expected, identical to those depicted in Fig. 4. This is 
true not only for the oscillations which are periodic in 
T/2, but also for the global enveloping functions. The 
scaling with A'^/^ is also observed for the dipole accelera- 
tion (Fig. |l^(b)). The parameters used in the figure are 
typical for quantum wells and solid-state systems M . 
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FIG. 13. Harmonic spectrum for the same parameters as in 
Fig. 4 {X — 1), compared to those obtained for several field 
strengths Eo, field frequencies u, transition frequencies ojio 
and matrix dipole elements xio, chosen such that 71 = 0.0391 
and 72 — 0.3197, i.e., the same as in Fig. 4. These parameters 
are displayed in Table IV. Part (a) shows the whole spectra, 
whereas part (b) displays both spectra for harmonic order 
10 < Af < 20, such that their substructure can be seen. The 
field is switched on linearly within two cycles. 
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Another interesting aspect concerns the resulting har- 
monic spectra. Even though, in absolute terms, these 
spectra have different cutoff frequencies and different 
global intensities, for equal 71 and 72 they have the same 
shape. Not only the number of harmonics is the same. 
In addition, all substructure in the spectra looks strik- 
ingly similar. These features can be easily understood: 
the global intensity decrease is related to the decrease 
in amplitude of the dipole acceleration and the identical 
shapes are a consequence of the fact that the populations 
of the adiabatic states, as well as all oscillations present in 
the dipole acceleration, remain invariant under the scale 
transformations discussed here. This is shown in Figs. 
p^(a) and (b), for several dilatation factors A. The corre- 
sponding field and two-level atom parameters are given 
in Table 0. 
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FIG. 14. Global structures as functions of time, for; (a) the 
populations |C;f of the adiabatic states; (b) the dipole 
acceleration x{t). The field strength, the field frequency, the 
transition frequency and the dipole matrix element were cho- 
sen as Eo = 0.62 a.u., lu — 0.05 a.u., loio = 0.409 a.u. and 
xio = 1.066 a.u., respectively. These parameters are slightly 
diflerent from the ones in Fig. 4, but give 71 = 0.0378, 
72 = 0.3094. The field is switched on linearly within two 
cycles. The time is given in units of the field cycle. 
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FIG. 15. Harmonic spectrum for the same parameters as 
in Fig. 4, compared to the one obtained for Eo = 0.62 a.u., 
w = 0.05 a.u., LOio = 0.409 a.u. and xio — 1.066 a.u., re- 
spectively. These parameters give 71 — 0.0378, 72 ~ 0.3094, 
whereas the ones in Fig. 4 yield 71 = 0.0391, 72 ~ 0.3197. 
Part (a) shows the whole spectra, whereas part (b) displays 
both spectra for harmonic order 10 < A'^ < 20, such that their 
substructure can be seen. The field is switched on linearly 
within two cycles. 

On the other hand, the behavior of the system can al- 
ready be altered by small variations in 71 and 72. For 
instance, in Fig. |lj we consider a slightly larger field am- 
plitude than in Fig. 4, which give different 71 and 72. In 
this case, one observes a radically different pattern for the 
populations \C^{t)\'^ [Fig. [l^(a)] and the dipole acceler- 
ation [Fig. |lj(b)]. As a direct consequence, the spectra 
does not exhibit the same substructure [Figs. |l^(a) and 
(b)]. 



VI. CONCLUSIONS 

The results discussed in the previous sections lead 
to the main conclusion that the three-step model and 
the two-level atom are not completely different physi- 
cal pictures for describing high-harmonic generation, as 
commonly believed. Indeed, in both models, this phe- 
nomenon takes place as a result of a three-step process. 
Hints that a correspondence between both physical pic- 
tures might exist have been provided in the literature 
We go however beyond such studies, giving evi- 
dence that a three-step mechanism exists in the two-level 
atom case and analysing its features in detail. 

In the usual form of the three-step model, there is pop- 
ulation transfer from the atomic ground state to a state in 
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the continuum, i.e., tunneling or multiphoton ionization. 
The electron then propagates in the continuum within a 
time interval t — ti — gaining a certain amount of 
kinetic energy which is converted into harmonic radia- 
tion at a time ti , when there is population transfer from 
the continuum to the ground state, i.e., recombination. 
In the two-level atom framework, a very similar process 
takes place: there is population transfer from the field- 
dressed state \(f>o{t)) to the state \4>ti't)) ^ time to 
for which an avoided crossing occurs. Subsequently, the 
system acquires energy from the field within the interval 
T = ti — to, and, at a further time ti, when population 
transfer from back to |0c?(i)) takes place, this 

energy is released in form of harmonic radiation. Thus, 
the main difference between the three-step model and the 
two-level atom physical pictures is that in the latter case, 
the three steps do not involve a continuum state, but a 
field-dressed bound state. 

Further similarities are observed in the time profile of 
high-harmonic generation. In both cases, the population 
transfers which contribute to the generation of a partic- 
ular set of harmonics occur at very specific times. In the 
usual three-step model, these times are such that the en- 
ergy of a particular harmonic must be equal to the sum 
of the kinetic energy of the electron upon return and the 
atomic ionization potential. The same line of argumen- 
tation holds in the two-level case, but now the harmonic 
energy must be equal to the energy difference between 
the adiabatic states at these times. 

Specifically for monochromatic driving fields, both 
models share several features. Both in the three-step 
model and in the two-level atom case, there is a single 
time corresponding to the generation of the cutoff har- 
monic. In the former model, this time corresponds to 
the maximal kinetic energy the electron may have, upon 
return, whereas in the latter model, it gives the maxi- 
mal energy difference between the adiabatic states. Also 
for both cases, this time splits into two sets of times as 
the harmonic energy decreases. The constructive inter- 
ference between the corresponding population transfers 
originates the plateau in the high-harmonic spectra. This 
pattern repeats itself every half cycle of the driving field. 
This is a direct consequence of the periodicity of the rel- 
evant physical quantities, namely the electron kinetic en- 
ergy in the three-step model and the adiabatic ener- 
gies in the two-level atom case. All these features are 
observed as peaks in the Gabor transform of the dipole 
acceleration. In the three-step model framework, analo- 
gous studies have been performed in 

Also for bichromatic driving fields, there are several 
characteristics which are present in both models. A good 
example is the multiple cutoff structure. Indeed, the har- 
monic spectra in this case may exhibit several cutoffs, 
which, depending on the model in question, are given by 
the maxima of either the electron kinetic energy or of 
the energy difference between the adiabatic states. The 
number of these cutoffs, as well as their energies or the 
corresponding population-transfer times, are determined 



by the frequency ratio n, the field-strength ratio C and 
the relative phase 9. For both the three-step model and 
the two-level atom, all peaks in the Gabor spectra can 
be traced back to the population-transfer times. In one 
or the other case, these population transfers occur either 
between the adiabatic states (Sec. IV), or between the 
ground-state and the continuum 

Similarities between the two models are also observed 
for the probability that the "first step", i.e., population 
transfer, takes place. In the three-step model, this prob- 
ability, per unit time, is roughly given by the quasi-static 
tunneling rate V - exp[-C/\E{to)\] |3|. A strong field 
E{to) at the ionization time to yields strong harmonics at 
the recombination time ti . This relation is very useful for 
controlling harmonic spectra, as for instance the relative 
intensities of a double plateau (see, e.g., ||,^ for concrete 
examples). Within the two- level atom framework and in 
the monochromatic case, to first approximation, the field- 
dependent terms of the two-level Hamiltonian can be lin- 
earized at the crossings Q. Thus, the population trans- 
fer between the exchanged states can be computed by 
means of the Landau-Zener model (2^j3l[] . This probabil- 
ity is approximately given by 7-" ~ exp[— C'7r/(2a::io£'o)], 
such that the Rabi frequency, in the two-level atom, plays 
a similar role as E{tQ) in the three-step model. In gen- 
eral, however, there is not always a simple expression for 
the population transfer at a level crossing j2^j3^, such 
that V has to be computed according to the problem at 
hand. For instance, V may be rather complicated for 
bichromatic fields. This is a limitation for controlling 
high-harmonic spectra in this latter case. 

A particularity of the two-level atom is that the very 
same distortions caused by the additional field in the 
field-dressed energies, as functions of time, are also 
present in the adiabatic-state populations |C^(t)P and 
in the dipole acceleration. Specifically for the bichro- 
matic field addressed in this paper, i.e., a, lo — 2uj field, 
the whole pattern is no longer periodic in T/2, but in T. 
This is a consequence of the periodicity of the adiabatic 
states, which changes with the additional driving wave. 
A similar feature occurs in the three-step model frame- 
work, due to an analogous change in the electron kinetic 
energy upon return (see, e.g., j^J^ for a discussion of this 
issue). 

An interesting issue which is not discussed in this paper 
concerns the influence of ionization or feedback mecha- 
nisms on the time profiles of harmonic generation by a 
two-level atom. In a previous paper it was shown that 
the main contributions to harmonic generation from a 
two-level atom whose states decayed according to quasi- 
static ionization rates occurred at minimal field. These 
results did not agree with the bound-bound transitions 
computed from the numerical solution of the Schrodinger 
equation for a gaussian potential with two strongly cou- 
pled bound states . The strikingly different time pro- 
files obtained in the present paper for HHG in a closed 
two-level atom suggest, however, that these features are 
stongly influenced by ionization. Therefore, more accu- 
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rate descriptions of ionization and an adequate feedback 
mechanism from the continuum would be necessary in 
the two- level atom case with unstable levels. The influ- 
ence of level widths on the population transfer between 
quantum states is discussed in | |36|] . 

Finally, there are scaling laws which allow extending 
the studies performed in this paper to a broader param- 
eter range. In fact, we have shown that the important 
parameters for determining the physical behavior of the 
system are 71 = u}/{2xwEo), and 72 = wio/(2a;ioSo), 
which denote the ratio of the field and transition fre- 
quencies to the Rabi frequency, respectively. As long as 
71 and 72 are kept constant, driving fields of completely 
different strengths and frequencies acting on systems of 
completely different energy gaps can yield similar spec- 
tra. For bichromatic fields, an additional requirement 
for this invariance are fixed field-strength ratio field- 
frequency ratio n and relative phase 0. 

A concrete example of a system for which these prop- 
erties may be applied is for instance a quantum well with 
LOiQ ~ 10"'' a.u., and xio ~ 100 a.u., subject to a field 
of strength £^0 ~ 10"^ a.u. and frequency cj ^ 10~^ 
a.u. [ p4| . Transitions between two subbands in these sys- 
tems are described very frequently by the semiconduc- 
tor Bloch equations in the Hartree-Fock approximation 
psi . In case collective effects can be neglected, the corre- 
sponding Hamiltonian reduces to a two-level one-particle 
Hamiltonian. In such a case, the results of the present 
paper are expected to be applicable. 
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TABLE I. Level-crossing times to, population transfer 
times ti and the corresponding harmonic energy f2, for the 
parameters of Fig. 2. The times are given in units of the pe- 
riod T = 2-K juj. The harmonic orders, together with the ap- 
proximate harmonic energies in units of the cutoff frequency 
Am, are given in the remaining two columns. This pattern 
repeats itself every half-cycle of the driving field. 





ti/T 


harmonic order 




0.5 


0.25 


43 


1 


0.5 


0.14 0.36 


35 


0.8 


0.5 


0.09 0.41 


25 


0.6 


0.5 


0.05 0.45 


17 


0.4 
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TABLE II. Times for the population transfers between the extrema of the adiabatic states, with the approximate order of 
the corresponding cutoff harmonic, for a bichromatic field given by Eq. ([l7|), with relative phase 9 = and several field-srength 
ratios = Eo2/Eoi. The field and two-level atom parameters are the same as those used in Fig ^. No entry means that 
corresponding maxima do not exist. This pattern repeats itself every cycle T = 2n/u! of the driving field. 



to/T 


C = 0.2 

tui/T 




to/T 


C = 0.5 

tlM/T 




to/T 


C = 0.8 

tlM /T 







0.20 


43 





0.17 


43 





0.15 


43 


0.5 


0.80 


43 


0.5 


0.83 


43 


0.36 


0.42 
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0.85 


43 














0.64 


0.58 


9 



TABLE III. Times for the population transfers between the extrema of the adiabatic states, with the approximate order 
of the corresponding cutoff harmonic, for a bichromatic field given by Eq. dl^), with relative phase 6 = n/2 and several 
field-srength ratios ^ = -E02/-E01. The field and two-level atom parameters are the same as those used in Fig No entry 
means that corresponding maxima do not exist. This pattern repeats itself every cycle T = 2n/Lj of the driving field. For 

= 0.8, there are additional avoided crossings at 0.25rmod T. 
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tlAl/T 


^Im/c^ 
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tlM/T 
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tlM/T 
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0.53 
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43 
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1.25 


30 
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1.08 


23 
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1.05 


24 








0.94 
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1.45 


24 
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1.45 


24 
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TABLE IV. Field and two-level atom parameters, given in 
atomic units, together with the dilatation factor A. All param- 
eters have been chosen such that 71 = 0.0391 and 72 = 0.3197. 





Eq 
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OJlO 


A 


1.066 


0.6 


0.05 


0.409 


1 


9.535 
47.673 


8.385 X 10~* 
6.71 X 10"^ 


6.25 X 10"" 
2.5 X lO"'^ 


5.1125 X 10"-' 
2.045 X 10"* 


1/80 
1/2000 



